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We report on a search for new resonant states in the process 77 —> DD. A candidate C-even 
charmonium state is observed in the vicinity of 3.93 GeV je?. The production rate and the angular 
distribution in the 77 center-of-mass frame suggest that this state is the previously unobserved Xc 2 , 
the 2 ®P 2 charmonium state. 

PACS numbers: 13.66.Be, 14.40.Gx 


The masses and other properties of the ground and 
excited states of charmonium provide valuable input to 
QCD models that describe heavy quarkonium systems. 
To date, radial excitation states of charmonium are es¬ 
tablished only for the = ^Si {ip) and, recently, 

the ^So {rjc) states. Although the lowest ^Pj states 
iXej) are already well established, their radial excitations 
have not yet been observed. 

The first radially excited Xcj states are predicted to 
have masses between 3.9 and 4.0 GeV/c^ y, which 
is considerably above DD threshold. If the masses of 
these states lie between the DD and D*D* thresholds, 
the Xco( 2 P) (Xco) and Xc 2 ( 2 P) (Xc 2 ) are expected to de¬ 
cay primarily into DD, although the Xc 2 could also decay 
to DD* if it is energetically allowed. (The inclusion of 
charge-conjugate reactions is implied throughout this pa¬ 
per.) Recently, two new charmonium-like states in this 
mass region, the A(3940) 0 and y(3940) 0, were re¬ 
ported by Belle. Neither of these states has been observed 
to decay to DD 0. 

In this paper we report on a search for the x'cj {J = ^ 
or 2 ) states and other C-even charmonium states in the 
mass range of 3.73 - 4.3 GeV/c^ produced via the process 
77 —> DD. 

The analysis uses data recorded in the Belle detec¬ 
tor at the KEKB e+e“ asymmetric-energy (3.5 on 8 
GeV) collider 0. The data sample corresponds to an 
integrated luminosity of 395 fb“^, accumulated on the 
T(45') resonance (^/s = 10.58 GeV) and 60 MeV be¬ 
low the resonance. We study the two-photon process 
e+e” ^ e~^e~DD in the “zero-tag” mode, where neither 
the final-state electron nor positron is detected, and the 
DD system has very small transverse momentum. 


A comprehensive description of the Belle detector is 
given elsewhere 0 . Charged tracks are reconstructed 
in a central drift chamber (GDC) located in a uniform 
1.5 T solenoidal magnetic field. The 2 : axis of the detec¬ 
tor and the solenoid are along the positron beam, with 
the positrons moving in the —z direction. Track trajec¬ 
tory coordinates near the collision point are measured by 
a silicon vertex detector (SVD). Photon detection and 
energy measurements are provided by a CsI(Tl) elec¬ 
tromagnetic calorimeter (ECL). Silica-aerogel Cherenkov 
counters (ACC) provide separation between kaons and 
pions for momenta above 1.2 GeV/c. The time-of-flight 
counter (TOE) system consists of a barrel of 128 plastic 
scintillation counters, and is effective for K/tt separation 
for tracks with momenta below 1.2 GeV/c. Low energy 
kaons are also identified by specihe ionization {dE/dx) 
measurements in the GDC. 

Kaon candidates are separated from pions based on 
normalized kaon and pion likelihood functions obtained 
from the particle identiheation system {Lk and re¬ 
spectively) with a criterion, Lk/{Lk + > 0 . 8 , which 

gives a typical identification efficiency of 90% with a 
probability of 3% for a pion to be misidentihed as a kaon. 
All tracks that are not identified as kaons are treated as 
pions. 

Signal candidates are triggered by a variety of track- 
triggers that require two or more GDC tracks with as¬ 
sociated TOE hits, ECL clusters or a minimum sum of 
energy in the ECL. For the four and six charged track 
topologies used in this analysis, the trigger conditions 
are complementary to each other and, in combination, 
provide a high trigger efficiency, (96 ± 3)%. 

We search for exclusive DD production in the following 
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The symbols in parentheses are used to designate each 
of the final states. For the four-prong processes (N4 and 
N5) the selection criteria are: four charged tracks, each 
one with (L) a transverse momentum to the z axis in the 
laboratory frame of pt > 0.1 GeV/c; two or more tracks 
must have (S) pt > 0.4 GeV/c and 17° < 6 < 150°, 
where 9 is the laboratory frame polar angle; no pho¬ 
ton clusters with an energy greater than 400 MeV; the 
charged track system consists of a iF+iF“7r+7r“ combi¬ 
nation; the K^tt^ combination with the larger invariant 
mass should lie within ±15 MeV/c^ of the nominal 
mass. For the N4 process, we require that the 
combination with the smaller invariant mass be within 
^20 MeV !(? of the nominal mass. For the N5 process, 

we require that the remaining Kt: combination, when 
combined with a 7r° candidate, has an invariant mass in 
the range 1.83 GeV/c^ < M(iF+7r“7r°) < 1.89 GeV/c^. 
Candidate 7r*^’s are formed from pairs of photons with 
energies greater than 20 MeV, which fit to the 7r° ^ 77 
hypothesis with < 4. If there are multiple 7r° can¬ 
didates, we select the one that results in M(iF+7r“7r°) 
closest to the nominal mass. 

For the six-prong processes (N6 and C6), we re¬ 
quire exactly six tracks with particle assignments 
iF+iF“7r+7r“7r+7r“, where all six pass the looser track 
criteria, indicated by (L) above. In addition, either two 
to four tracks must pass the more stringent track criteria 
(S) or at least one track has pt > 0.5 GeV/c and the 
sum of ECL cluster energies is less than 0.18-\/s, where 
the cluster energies are measured in the e^e“ center-of- 
mass (c.m.) system. For the N6 process, one combina¬ 
tion is required to have |AM|i = — m£)o| < 

15 MeV/c^ while the remaining tracks have |AM |2 = 
|M(Ar“7r+7r“7r+) — TO£|o| < 30 MeV/c^. When there 
are multiple combinations, we choose the one with the 
smallest |AM|i ± |AM| 2 . For the C6 process, we re¬ 
quire — m£)+| < 30 MeV/c^ for each of 

the charge combinations, where m£)+ is the nominal D'^ 
mass. 

For all processes, we require that there be no ex¬ 
tra 7r° candidates with transverse momenta larger than 
100 MeV/c. We also apply the following kinematical 
requirement to the DD candidate system: Pz{DD) > 
{M{DDf -_49 GeVVc'‘)/(14 GeV/c^) ± 0.6 GeV/c, 
where Pz{DD) and M{DD) are the momentum com¬ 
ponent in the z direction in the laboratory frame and 
the invariant mass, respectively. This condition removes 
events from initial-state radiation (ISR) processes, such 
as e+e" —> in which the photon is emit- 
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FIG. 1: Invariant mass distributions of (a) in N4 

candidate events, (b) in N5 candidate events, (c) 

in N 6 candidate events, (d) in N 6 can¬ 
didate events, and (e) in C 6 candidate events. An 

accompanying D meson candidate is required in each event 
sample. 

ted in the forward direction with respect to the incident 
electron. We compute M{DD) using the measured 3- 
momenta of each D candidate (Pd) and energy deter¬ 
mined from Ejj = P^ + m?j ^, where mo the nominal 
mass of the neutral or charged D. 

The invariant mass distributions for D meson can¬ 
didates reconstructed with the above requirements are 
shown in Fig. 1. 

We calculate Pt{DD), the total transverse momen¬ 
tum in the e^e“ c.m. frame with respect to the inci¬ 
dent e“'"e“ axis that approximates the direction of the 
two-photon collision axis. We apply the requirement 
Pt{DD) < 0.05 GeV/c to enhance exclusive two-photon 
77 —> DD production. In the invariant mass region 
M{DD) < 4.3 GeV/c^, we find 86 N4-process events, 
60 N5-process events, 168 N 6 -process events and 128 C 6 - 
process events. 

In Figs. 2(a) and 2(b) we show the M{DD) distribu¬ 
tions separately for (sum of N4, N5 and N 6 ) and 

D~^D~. The invariant-mass distribution for the com¬ 
bined D^D^ and D^D~ channels is shown in Fig. 2(c). 
There, two event concentrations are evident: one near 
3.80 GeV !(? rather close to the threshold of DD and an¬ 
other near 3.93 GeV/c^. Each distribution of the four 
decay combination modes shows an enhancement near 
the latter invariant mass. We apply an unbinned maxi¬ 
mum likelihood fit to the combined data in the region 
3.80 GeV/c^ < M{DD) < 4.20 GeV/c^ using a rel¬ 
ativistic Breit-Wigner signal function for the resonant 
peak near 3.93 GeV/c^ plus a background of the form 
M{DD)~°‘, where a is a free parameter. The invari¬ 
ant mass dependence of the efficiency (decreasing by 
10% for an increase of the invariant mass from 3.80 to 
4.20 GeV/c^) and the two-photon luminosity function 
are taken into account in the resonance function. These 
are computed using the TREPS Monte-Carlo(MC) pro¬ 
gram for e“'"e“ ^ e'^e~DD production together with 
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JETSET7.3 decay routines 0 for the D meson decays 
(using PDG2004 [ij values for the decay branching frac¬ 
tions). We find from the MC study that the product of 
the efficiency and branching fractions of the two D decay 
modes in the D'^D~ channel is about 50% of that in the 
channel. 

The results of the fit for the resonance mass, width 
and total yield of the resonance are M = 3929 ± 
5{stat) MeV/c^, E = 29 ± 10{stat) MeV and 64 ± 
18{stat) events, respectively. The mass resolution, which 
is estimated by MC to be 3 MeV/c^ is taken into account 
in the fit. The statistical significance of the peak is 5.3a, 
which is derived from -^2 ln(Tinax/.^o)j where Lmax and 
Lq are the logarithmic-likelihoods for fits with and with¬ 
out a resonance peak component, shown in Fig. 2(c) as 
solid and dashed curves, respectively. 

Systematic errors for the parameters M and E are 
2 MeV/c^ and 2 MeV, respectively. The former is par¬ 
tially due to the uncertainties on the D meson masses 
(1 MeV/c^ for the resonance mass). We also consider 
the effect of choosing different Breit-Wigner functional 
forms for spin 0 and 2 resonances and wave functions in 
this error. 

The Pt{DD) distribution in the peak region, 
3.91 GeV/c^ < M{DD) < 3.95 GeV/c^, is shown in 
Fig. 3. Here the Pt requirement has been relaxed. The 
experimental data are fitted by a shape that is expected 
for exclusive two-photon DD production plus a linear 
background. We expect non-charm and non-exclusive 
backgrounds to be nearly linear in Pt{DD). The fit 
uses a binned-maximum likelihood method with the nor¬ 
malizations of the two components treated as free pa¬ 
rameters. The linear-background component, 1.8 ± 0.6 
events for Pt{DD) < 0.05 GeV/c^, and the goodness 
of fit, x^/d.o.f = 14.2/18, indicate that the events in 
the peak region originate primarily from exclusive two- 
photon events. 

The Pt{DD) distribution produced by DD* and D*D* 
events is expected to be distorted by the transverse mo¬ 
mentum of the undetected slow pion(s), which peaks 
around 0.05 GeV/c (dashed histogram in Fig. 3). Such a 
distortion is not seen in the observed Pt distribution. 

We investigate possible backgrounds from non-DD 
sources using H-sideband events. The histogram in 
Fig. 2(c) shows the invariant mass distribution for events 
where the H-meson is replaced by a hadron system from a 
D-signal mass sideband regions above and below the sig¬ 
nal region with the same width as the signal mass region. 
Here we use two types of sideband events: one where one 
H-meson candidate is in the signal mass region, and an¬ 
other where both entries are from the sidebands. Since 
there is no significant event excess in the former type over 
the latter, we conclude that the sideband events are domi¬ 
nated by non-charm backgrounds. We combine them and 
appropriately scale in order to compare to the DD signal 
yield. We conclude that the candidate events are domi- 



FIG. 2: Invariant mass distributions for the (a) chan¬ 

nels and (b) the D'^D~ mode, (c) The combined M(DD) 
distribution. The curves show the fits with (solid) and with¬ 
out (dashed) a resonance component. The histogram shows 
the distribution of the events from the D-mass sidebands (see 
the text). 


nated by DD (inclusive or exclusive) events in the entire 
mass region. 

Figures 4(a) and (b) show the M{DD) distributions 
for events with |cos0*| < 0.5 and |cos 6 **| > 0.5, respec¬ 
tively, where 9* is the angle of a Z1 meson relative to 
the beam axis in the 77 c.m. frame. It is apparent that 
the events in the 3.93 GeV/c^ peak tend to concentrate 
at small |cos0*| values. The points with error bars in 
Fig. 4(c) show the event yields in the 3.91 GeV/c^ to 
3.95 GeV/c^ region versus |cos0*|. Background, esti¬ 
mated from events in the M{DD) sideband, is indicated 
by the histogram. The solid curve in Fig. 4(c) shows 
the expectation using sin^ 9* to represent the signal from 
a spin -2 meson produced with helicity -2 along the inci¬ 
dent axis . A term proportional to 1 -I- a cos^ 9* 

that interpolates the background (dotted curve) is also 
included. A small nonuniformity of the signal acceptance 
in the c.m. angle is taken into account. The comparison 
to the data has y^/d.o.f. = 1.9/9. Here the functions 
are normalized to the numbers of signal and background 
events obtained from the fit of the invariant mass dis- 
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Pt(DD) (GeV/c) 


FIG. 3: Experimental Pt{DD) distribution (points with error 
bars) for events in the 3.91 < M(DD) < 3.95 GeV/c^ region 
and the fit (histogram) based on the exclusive 77 —> DD 
process MG plus a linear background (dotted line). The dot- 
dashed line shows the location of the Pt selection requirement. 
The dashed histogram shows the expected distribution of the 
77 —> DD* process followed by D* —>■ Du with an arbitrary 
normalization. 



FIG. 4: M{DD) distributions for (a) |cos0*| < 0.5 and 
(b) I cos 0*1 > 0.5. (c) The |cos0*| distributions in the 

3.91 < M{DD) < 3.95 GeV/c^ region (points with error bars) 
and background scaled from the M{DD) sideband (solid his¬ 
togram). The solid and dashed curves are expected distribu¬ 
tions for the spin two (helicity two) and spin zero hypotheses, 
respectively, and contain the non-peak background also shown 
separately by the dotted curve. 


tribution, 46 and 33 events, respectively. A compari¬ 
son using a constant term to represent the signal from 
a spin -0 meson (dashed curve) gives a much poorer fit: 
X^/d.o.f. = 23.4/9. The data significantly favor a spin 
two assignment over spin zero. 

No charmonium state that decays into DD with a mass 
near 3.93 GeV/c^ has been previously reported. This 
observation cannot be attributed to a new = 1 
meson (ip) produced by ISR processes, because there are 
no structures as large as the signal at this mass in the 
e“'"e“ hadronic cross section. 

Using the number of observed signal events, the 
branching fractions and efficiencies for the four de¬ 
cay channels, we determine the product of the two- 
photon decay width and DD branching fraction to be 
r.^.^(Z(3930))S(Z(3930) ^ DD) = 0.18 ± 0.05(stat) ± 
0.03(s?/s) keV, assuming production of a spin-2 me¬ 
son. Here, we define U(Z(3930) ^ DD) = 

S(Z(3930) ^ D°D°) -be(Z(3930) ^ D+D-) and as¬ 


sume 6(Z(3930) ^ D+D-) = 0.89S(Z(3930) ^ D°D°) 
according to isospin invariance and including the effect of 
the mass difference between D^ and D^ mesons, where 
Z(3930) is used as a tentative designation for the ob¬ 
served state. 

We assign a 17% total systematic error to the mea¬ 
surement of the product of the two-photon decay width 
and the branching fraction, as shown in the above result. 
This is primarily due to uncertainties in the track recon¬ 
struction efficiency (7%), selection efficiency ( 8 %), kaon 
identification (4%), choice of the fit function and back¬ 
ground shape (5%), luminosity function (5%), and the 
U-meson branching fractions (9%), added in quadrature 
with other smaller factors. 

The observed signals for the D^D^ and D^D~ modes 
are consistent with isospin invariance. The ratio of the 
branching fractions is measured to be ,B(Z(3930) ^ 
D+D-)/B(Z(3930) D°D°) = 0.74 ± 0A3(stat) ± 

0.16(sys). The results on mass, decay angular distribu¬ 
tions and T-yjB(—>- DD) are all consistent with expecta- 
tions for the x'c 2 ^ ^^e 2^P2 charmonium state |2Ll3lll3l|. 

In summary, we have observed an enhancement in DD 
invariant mass near 3.93 GeV/c^ in 77 —> DD events. 
The statistical significance of the signal is 5.3tT. The ob¬ 
served angular distribution is consistent with two-photon 
production of a tensor meson. Results for the mass, 
width, and the product of the two-photon decay width 
times the branching fraction to DD are: M = 3929 ± 
3(stat)±2(sys) MeV/c^, T = 29±10(stat)±2(sj/s) MeV 
and T^^B(-^ DD) = 0.18 ± 0.05(stat) ± 0.03(sj/s) keV 
(assuming J = 2), respectively. The measured proper¬ 
ties are consistent with expectations for the previously 
unseen -x!c 2 charmonium state. 
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